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I. INTRODUCTION 



All-optical nonlinearities at low light levels are instrumental to the optical implementation of quantum information 
processing systems. Large nonlinearities can be achieved when the light-matter interaction is resonant, in which case 
however absorption typically plays a detrimental role. A great deal of attention has thus been devoted to a variety of 
schemes based on the regime of electromagnetically induced transparency (EIT) [l[ in which this counterbalance can 
be overcome. In a typical EIT configuration, a strong control field, coupling two unpopulated levels of a A system, 
creates a transparency window for a weak probe beam. This basic scheme was extended by including additional atomic 
levels coupled by laser or microwave fields in the double- A 0|, tripod 0,0] or inverted-F Q configurations, 

just to name a few. Under EIT-like conditions the weak interaction between two photons (or weak classical pulses) may 
become largely enhanced. In particular, very large cross nonlinear effects have been predicted leading to new types of 
polarization phase gates in an optically dressed medium in the M 0], tripod 10 . 11 | or inverted-F 12 1 configuration. 



Relevant experimental work has very recently been reported in Ref. [13|. In these instances, a significant change 
of the phase of one of the propagating pulses is achieved due to the cross-Kerr effect induced by the other pulse. 



Analogous effects have recently been studied, both theoretically and experimentally, for an inverted-F system IJ], 
as well as for a four- level A^-type [l3| and five- level system 16 1. Important developments concern also the dynamic 
control of the process leading, e.g., to light slowdown, storage and release 17|, |l8|. These results are expected to pave 
the way towards constructing all-optical logical devices. 

Unlike most typical EIT configurations employing a running wave coupling field, using a control standing wave 
coupling beam opens the possibility of creating an ail-optically tunable Bragg mirror. In such a novel kind of 
dynamically tunable metamaterial 1^ IB] an incoming probe light encounters a spatially periodic optical structure 
and is subject to Bragg scattering. For specific frequency ranges transmission and stop bands appear, the properties 
of which can be steered by a proper choice of the control field. The field-induced dispersion of the medium and the 



transmission and reflection spectra have been analyzed, e.g., in Refs. |19l - l2ll |. including the case of a quasi-standing 
wave coupling field 0. It was also possible to stop and store a pulse inside the medium, and then retrieve it in the 
form of a stationary light pulse which could not leave the medium due to the standing wave character of the releasing 
control beam 1^, 23 1. It has also been shown that high nonlinearities can be achieved for stored light pulses 24 1. An 
all-optical dynamic cavity for the confinement of light pulses based on standing wave EIT Bragg mirrors has recently 
been proposed [ist . 

Very recently, we have pointed out that in the presence of a trigger beam in a tripod configuration, a probe partially 
reflected from the periodic structure induced by a strong standing wave EIT coupling field undergoes large cross-Kerr 
nonlinear phase shifts. Such a novel configuration is specifically apt for developing a phase-tunable beam splitter in 
which both the amplitudes and the phases of the trasmitted and reflected beam can be controlled. We have given some 
illustrative results valid in a narrow range of frequencies owing to the restrictive assumption that the atomic population 
is symmetrically distributed between the two lower levels that are coupled with the upper level respectively by the 
probe and trigger pulses (26l | . Here, we provide a comprehensive theoretical study of the problem, including a proper 
treatment of the population redistribution. The latter point turns out to be crucial to extend the range of possible 
probe detunings and, thus, optimize the control possibilities. Our approach combines analytical methodologies used to 
describe propagation effects in both nonlinear and spatially periodic media, and allow us to numerically examine the 
phase shifts of the reflected and transmitted probe fleld induced by the trigger's presence over a wide range of probe 
detunings as well as their dependence on various parameters characterizing both the atomic medium and the three 
laser fields. Our results indicate that the tripod standing wave EIT configuration makes control over the cross-Kerr 
effect more versatile than for a running wave EIT configuration l^, 11 1 in which no reflected beam appears. 



The paper is organized as follows. In Section II, we present the basic theoretical approach to calculate the cross- 
Kerr effect in a medium additionally dressed by a quasi-standing wave EIT coupling fleld, based on expanding the 
medium susceptibility into a power series with respect to the probe and trigger fields and into a Fourier series 
appropriate to the spatial periodicity imposed by the coupling field. In Appendix A we derive in detail the formulae 
for the susceptibility's expansion coefficients, while in Appendix B we present the way to evaluate the population 
redistribution. In Section III, we present and discuss the numerical results for the trigger induced nonlinear phase 
shifts of the refiected and transmited probe fields as functions of the probe detuning for various combinations of 
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control field intensities, relaxation rates and sample length. Finally, we draw our conclusions. 



II. THEORY 

We consider a tripod system driven by a strong control field {£c) coupling the ground level 2 and the upper level 
0, and by two weaker laser fields called probe {£p) and trigger {£t), coupling with 1 and 3, respectively, as shown in 
Fig. [TJ We express the fields through their complex amplitudes ep.c.t slowly varying in time: 

Each of the couplings is detuned by Sj = Ej + hujj — Eq, j = 1,2,3. Here Ej stands for the energy of the state j and 
LUj denotes the frequency of the j-th field. 

All the fields propagate along the z axis. In general the control wave is allowed to have two counter-propagating 
components of arbitrary ratio. In the case of a running control wave, the component antiparallcl to the incoming 
probe and trigger fields is zero. In the case of a perfect standing wave, both components are equal. Thus, the control 
field provides a spatial lattice of period ^2 being the field's wave number. In the configuration proposed below 
k2 may be slightly different from the probe and trigger fields' wave vectors fci,3, but that can be corrected by tilting 



both components of the control field by an angle 9 with respect to the z-axis so that W2 cosO/c « fc^ 3 |22|. 

Such a configuration can be realized in a cold gas of ^''Rb atoms. The states 1 and 3 may correspond, e.g., to 
Zceman sublevels |5S'i,F = 2,m = { — 1,1} >, state 2 to |55i,F = l,m = 1 >, and the upper state to the level 
\5P3,F — l,m = >. The scheme could be used as a polarization phase gate, as was first proposed in Ref. [lol |. 
Each of the two pulses interacts with the medium when it has the proper circular polarization: right for the probe 
and left for the trigger. When both pulses are properly polarized the nonlinear interaction between them gives rise 
to a cross-phase modulation. 

The Bloch equations for the atom + field system in the rotating wave approximation read: 
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uj-uik)t^ and g is the density matrix of the atoms in the Schrodinger 


= Sj - Sk ~ 


~ i7jfc, 7jfc = 7fcj, J, fc = 1, 2, 3, m = 0, 1, 2, 3. The Rabi frequencies are defined 



where amm = gmm: cToj ^ goje 
picture, A^o = Sj+ijjo, Ajk 

by ^p,c,t = —i2MMllE2£^^ where doj is the electric dipolc moment matrix element. 

The simplified model of relaxations adopted in the above equations takes into account spontaneous emission from 
the upper state to a lower state fc, described by the relaxation rates jQk and 7^.^. with 70A: — 5(711 + 722 + 733), 
as well as deexcitation and decoherence between the lower levels jjk-, j ^ k — 1,2,3, due to atomic collisions. The 
relaxation rates for the coherences ajk, j 7^ fc, and those for the collision-induced population relaxations were taken 
equal for simphcity (cf. Ref. [lol|). 

The propagation equation for the positive frequency part of the probe field, slowly varying in time, i.e. when 
|^-| << aji|ep|, reads in the Fourier picture, with the Fourier frequency variable identified with the probe detuning 



fLp[z,5x) ^ -'^Xp{z,5i)tp{z,5i), (2) 
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with the medium susceptibility for the probe given by 

Xp = - lim — — , (3) 

t-s-oo fteo Sip 

where A'' is the number of atoms per unit volume. The time limit means that we need to find the steady state solutions 
to the Bloch equations. 

The assumption that the probe and trigger are weak with respect to the control field |f2pp, |57tp << [r^d^ allows 
us to expand the susceptibility into the Taylor series: 

Xv-^x'p^+xf^\%\''+Xph^t?- (4) 

The explicit form of the Taylor expansion coefficients can be found in Appendix A. The term xl^' corresponds to the 
linear susceptibility, while the third-order terms Xpp and Xpt represent the self- and cross-Kerr effect, respectively. 
The latter is of particular importance, as we are interested in the cross-phase modulation. 
We now write the control field as: 

fic = nte'^^-"- + n-e"^^-"-, (5) 

where $7^ of constant values are its forward and backward propagating parts. As the field is periodic in space, so are 
the medium optical properties. Therefore we expand the susceptibilities for both probe and trigger into the Fourier 
series: 

oo 



n— — oo 



The pulses pro pag ating in a periodic medium acquire their reflected components. In the two-mode approximation 



(for details see 20|, |22|) we write: 

^]p,^ = ^]+,e'^^^ + ^]-,e-''=^^ (7) 

where fJ^ j are slowly varying in space. 

We now insert Eqs. (|6l7p into Eq. ([2|) and drop terms rapidly oscillating in space. To write the propagation 
equations in the final form we make use of the relation fc2 = ^ and set W2 = Wi -I- Awi. For the probe we find: 

= -^(xi7+ + yf7;), (8) 
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with Sp t = l^lp + |rip jp. The trigger equations are obtained by interchanging the indices {p, 1} <-> {i, 3}. Explicit 
formulae for the susceptibilities Fourier components as well as their detailed deriviation can be found in Appendix A. 

As we will change the probe detuning in a wide range while keeping the trigger detuning 5^ constant, wc introduce 
an asymmetry in the populations txn and eras even for ~ Additionally, the symmetry may be spoiled by the 
fact that El (the energy of the state 1) is slightly smaller than E^, and in a cold medium transitions from 3 to 1 
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arc possible while those from 1 to 3 are not. In previous works (see Refs. 10, 26|) the probe and trigger fields were 
only allowed detunings such that w ^3. In that case the stationary population distribution was assumed to be 
cii = 0^33 = \- Such an approach is simple, but it is well justified only in the resonance region of frequencies. As 
we show in the section IIIIl for a wider range of pulse detunings it is necessary to make a better estimation of the 
population distribution. In Appendix B we present the way to evaluate the populations in the case of a quasi-standing 
coupling beam ^ fi^. as in Ref.j22j). 



III. RESULTS AND DISCUSSION 



In this section we present numerical results illustrating the cross-Kerr effect in our system and its dependence 
on the system parameters. In particular, we discuss how the cross-phase modulation is affected by the population 
redistribution. 

We have solved the propagation equations (|8l9p with the boundary conditions corresponding to both probe and 
trigger originally propagating towards positive z: 



^pA^ = i,'5l,3) = 0, 



(10) 



where r2op,t describes the amplitude of the incoming probe/trigger pulse and L is the sample length. We have used 
the susceptibilities' Fourier components as presented in Appendix A. The level populations have been evaluated as 
described in Appendix B, using the initial amplitudes of the probe and trigger fields as input data. Note that, strictly 
speaking, the populations in the standing wave case are also modulated in space. However, as discussed below, the 
latter effect should not play a very important role in the case of an quasi-standing wave. 

The equations are solved iteratively until self-consistency is achieved, alternately for the probe and trigger, with 
the susceptibilities calculated using the fields obtained in the preceding step. In the first iteration for the probe the 
initial constant value of the trigger has been used. The numerical task in each step is thus solving a linear equation. 
We find two independent solutions with one-point boundary conditions and combine them to obtain the solution with 
the proper two-point boundary conditions given by Eqs. (jlOp . 

The transmission and reflection coefficients for the probe and trigger beams are defined by: 



Tp,ti^l,3) 
Rp,t{Sl,3) 



^ttiL,Si,3) 



^Op,t 



(11) 



no 



p,t 



while the phases of the transmitted and reflected fields are given by: 

^+,((5i,3) - arg(n+,(L,Ji,3)), 
Vp,ti.^i,i) = arg(rj-t(0,(5i,3))- 



(12) 



Unless otherwise specified, we use the following set of input data which are of order of those of typical atomic 
systems: fJ+ = 4 MHz, = 2 MHz, Opo = ^ta = 0.67 MHz, 710,20,30 = 0.67 MHz, 712,32,13 = 6.67 x lO"'' MHz, 
711,22,33 = 0.44 MHz, 82 = 6.67 MHz, 63 = 1.002 x 6.67 MHz, L = 1.06 mm, dio = ^30 = 8 x lO'^o C m, iV = 1.3 x lO^^ 
cm""^. Our choosing a quasi-standing wave is connected with the fact that for a perfect one the steady state population 
in the node regions is trapped in the state 2 due to relaxation effects and the medium becomes transparent to both 
probe and trigger. Therefore, it is reasonable to only consider the cases in which the control field at the quasinodcs 
is still strong enough to pump the population from the level 2. 

We first calculate the corrected values of the populations which are shown in Fig. [2] in a narrow (plot (a)) and 
wide (plot (b)) ranges of the probe detunings. Note that the populations cth and 0^33 may significantly differ from i. 
They may vary rapidly within the transparency window and the maximum value of G33 may be close to unity which 
apparently occurs when the frequency of the probe suits the energy interval between Eq and the energy of one of 
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the lower states dressed by ftc and J7t. Between the maxima there is a deep minimum of 1733 due to the population 
trapping in the state 1. The populations of (Too and (T22 are negligible, which is consistent with our assumption. 

If there are no relaxations in the system, then the steady state solution for the no-trigger case is trivial: 1733 = 1, 
o'li = Too = o'22 = due to spontaneous emission from the upper level to all the lower levels. There is no mechanism 
yet pumping the population out from 3. The nonncgligible relaxations provide such a mechanism and the plot 
illustrating the (733 dependence on the probe detuning is of a similar shape as in the trigger-present case (not shown). 

The phase shifts for a running control field (f2+ = 4 MHz, il^ ~ 0, the other parameters unchanged) are shown 
in Fig. [SK- We compare the results for the transmitted probe beam in the cases of equal or corrected values of the 
populations an and (T33 and for the trigger pulse present or absent. The values of the obtained phase shifts (p^{Si) 
differ significantly for the two approximations concerning the populations, except for three points: when the calculated 
population distribution is almost exact: an ~ (733 « 0.5, (T22 = c^oo ~ and within the resonance region 61 ~ 62 ~ S3. 
The trigger-induced (Kerr) phase shifts (the differences between the values given by pairs of curves) in the two cases 
are also different and there are frequency intervals in which the Kerr shifts are significantly lar ger in the case in which 
the populations have been calculated according to Eqs. (|B.3p . In this way the results of Ref. [10| can be generalized 
for a wider range of the probe detunings. 

In Fig. [31d we make a similar comparison for a reflected beam in the case of a quasi - standing control field (57+ = 4.19 
MHz, fl^ = 2 MHz, the other parameters unchanged). Again taking into account the population redistribution (lower 
plot) caused a qualitative change of the behaviour of the phase shifts as well as of their differences (the Kerr shifts) 
compared with the case of equal populations. The results are in agreement at the same three points as before. 
However, at 61 62 ^ S3, due to a numerical instability in the algorithm used to find the solutions of Eq. (jB.3| . 
the calculated density matrix is no longer positive definite and our results in this very narrow spectral region are 
unphysical (see the sharp spikes appearing in Fig. (2). 

In Fig. |3]we show the cross-Kerr phase shifts Aip^{Si) due to the trigger's presence for both the transmitted and 
reflected probe beams defined as Aip^ = (/^^(trigger on)— (^^ (trigger off). One can see that the shifts may be of order 
of one radian and there is a frequency range in which they do not change rapidly. The narrow minima or maxima in 
Fig. m correspond to the situation in which the phases ip^ vary rapidly. The present results show how our system 
behaves as an all-optically controlled Kerr medium. In particular, it may serve as a tunable beam splitter for the 
probe beam in which the phases of the reflected and/or trasmitted beam are significantly affected by the trigger 
beam. For example, for detunings of Si ^ 5.5 MHz or of Si ^ 6.4 MHz at which both the transmission and reflection 
coefficients are large (see Fig. [6]), the trigger induced phase shift of the trasmitted and reflected probe beams are 
both appreciable. On the contrary, for a detuning of Si 6.25 MHz at which both the transmission and reflection 
coefficients are considerable (see Fig. [6]), while the trigger induced phase shift for the transmitted probe beam is 
significant, that of the reflected beam is negligible. This beam splitting action accompanied by a trigger induced 
phase shift of the reflected and/or transmitted probe components is a unique feature of the tripod standing wave KIT 
configuration with respect to others previously considered. Taking into account the redistribution of the population 
among the atomic levels is however crucial. If the redistribution were ignored the phase shifts could only be estimated 
within rather limited ranges of frequencies. In particular, while for Si ~ 5.5 MHz the approximation an = CT33 = 0.5 
is still reasonable as shown in Fig. [21 for Si ^ 6.4 MHz or Si ~ 8.7 MHz it fails. In the following, in order to illustrate 
the dependence of the Kerr effect on various parameters, we discuss several numerical results obtained over a large 
detuning range using the method here developed to calculate the population redistribution. 

We now check how the phase shifts of both the transmitted field and reflected one depend on 'how much standing' 
the control field is. We change the right-propagating part of the control field $7+ while keeping the left-moving part 
constant: = 2 MHz. All the other parameters remain unchanged. The results are shown in Fig. [5^ for the 
transmitted beam and Fig. [Sb for the reflected one. All the plots cross in the resonance region where the phase shift 
is small. The phase shift reaches a maximal value for some Si which depends on $7+. The more intense is the 
control field the flatter is the plot and the larger is the peak's shift towards higher frequencies, which provides a way of 
controlling the phases. We find the same effect of the peak moving to the right when decreasing the left-propagating 
part of the control field fl~ while keeping il^ constant (not shown). The more intense is the control field the smaller 
is the trigger-induced phase shift. Around the resonance the reflected phase shift ip~ shows an oscillatory dependence 
on the probe detuning. The width of the frequency range where the oscillations arc present grows with 0+ but the 
number of peaks remains constant - again the phase shift's plot becomes flatter when the control field is increased. 
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The curves are complicated now but roughly we can say that the trigger-induced phase shift decreases when the 
control field becomes stronger. The extrema of the Kerr phase shifts A(/9~ in the intervals in which ip~ vary rapidly 
become less prominent and the distance between them increases for growing control fields. In general a too strong 
control field is not advantegous for generating considerable trigger-induced phase shifts. Then the results gradually 
turn into those typical of the usual EIT case. 

As expected, the transmission coefficient (see Fig. |6^) grows with f2+ and so docs the width of the transparency 
window. Again the plots corresponding to the reflection arc more complex, as shown in Fig. [6}3. In general the 
reflection coefficient decreases when f]+ is increased, but there are some frequency ranges where the dependence is 
more complicated. Yet if both parts of the control field are increased simultaneously, so that is constant, then 
not only the transmission but also the reflection coefficient increases for most of the probe frequencies. This is due to 
the absorption cancellation by a strong control field. 

We have also checked how our results depend on the relaxation rates jij , j ^ i — 1,2,3, due to interatomic collisions. 
The main observation is that the trigger-induced phase shift decreases for growing jij. On the other hand, as the 
relaxation rates grow, the frequency range where the reflected field's phase dependence is oscillatory becomes wider 
and the number of oscillations increases. As expected, the transmission and reflection coefficients decrease in general 
when the relaxation rates grow, but there are some frequency ranges where the reflection coefficient is a nonmonotonic 
function of jij. Populations ctoo and (T22 are negligible except for very large values of the relaxations (7i_, ^ 10~^ 
MHz) when they are of the order of a few percent. For increasing -fij the phase shifts become less steep and an 
approximation consisting in adopting constant values of the populations cth and 0-33 may be better justified. 

Increasing the length of the sample leads to increasing the trigger-induced phase shift of the transmitted field. For 
the considered lengths of the sample the nonlinear phase shift increases for growing L. However, for an extended 
sample both pulses become partially absorbed and their mutual impact is smaller. Elongating the sample leads to 
no significant change in the amplitude of oscillations in the refiected field's phase but it does influence the reflection 
coefficient in a nonmonotonic way, which provides another way of controlling the cross-Kerr effect. Thus, in order to 
have a considerable phase shift the coupling field cannot be too strong and the relaxation rates between the lower 
tripod states should not be too large; also the length of the sample should be appropriately optimized. 

IV. CONCLUSIONS 

We have presented a comprehensive study of the cross Kerr effect in the propagation of two weak laser fields, a 
trigger and a probe, in a medium of four-level atoms in the tripod configuration, dressed by a third strong coupling 
field in the quasi-standing wave configuration. This has been done by taking into the proper account the population 
redistribution. Using the relevant terms of the Fourier expansion of the linear and nonlinear susceptibilities, we have 
numerically solved the propagation equations for the right- (incoming) and left- (reflected) running components of the 
two weak flelds. In particular, the phase shifts of the transmitted and reflected probe induced by the trigger turn out 
to be as large as one radian. We have further shown that the population redistribution significantly affects the phase 
shifts of the transmitted and reflected beams unless very specific values of probe detuning arc used. Over a wide range 
of probe dctunings, on the contrary, we find that the trigger induced phase shifts for the reflected and transmitted 
probe can be flexibly and independently changed. We have flnally characterized the Kerr effect's dependence on the 
strength of the coupling fleld, the relaxation rates and the sample's length. 

In summary, by using a quasi-standing wave conflguration the amplitudes of the transmitted and reflected parts 
of the probe beam and their respective trigger-induced phase shifts can be controlled and optimized by acting on 
various parameters, and especially on the probe detuning. While previous work has shown the usefulness of the 
tripod conflguration in a running wave EIT regime to achieve a large cross-Kerr effect on the transmitted probe 
beam, our results extend such findings to the standing wave EIT regime in which both the reflected and transmitted 
probe beams can be separately controlled. 
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APPENDIX A 



The medium susceptibilities can be written as: 



iV|dioP CToi 
Xp = - hm — — , 

Xt = - hm — — . 

t->oo neQ lit 

and can be calculated with the use of the steady state solutions of the Bloch equations. 
From the last three of Eqs. ([T]) we find the steady state spin coherences: 

CT12 = , 

^12 

CT32 - , 

^32 

_ f^tCTio - f2*Cro3 

^13 = • 

^13 

We substitute these expressions into the stationary form of the equations 5 — 7 in Eqs. ([I} and thus obtain a set of 
three coupled equations for ctio, (J20, o'3o- As the control field is much stronger than any other field present in the 
system, it prevents the states and 2 from being populated. Therefore we set ctoo = tT22 = 0, as was done in Ref. [lo| . 
However, contrary to Ref. [l^l, here we do not assume the probe and trigger detunings to be almost equal. Hence the 
population distribution between the states 1 and 3 may be asymmetric and in general (Th 7^ (733. This crucial issue 
has been discussed in Section [HII Next we eliminate (T20 and arrive at: 



cr3o + —— T cro3 = "pCTn 



Vi3 L»Ai2A32 



-fToi = f^tO'33 



A31 LIA12A32 

A12 A32 

We combine the above expressions to obtain ^ and The solutions of the above set of equations expanded 



where D ~ A20 
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into the Taylor series have the form of Eq. (|4]), with: 

iVldiol^aii 



hen A10A12 - |Slc 



Xpn t, A * /" A A 



fco A*o(AioAi2-|f},|2)2' 



(3).. X _ ^l^iol 



A 



12 



1 



CTllA 



Ai3 A10A12 - If^cp 
12 , Cr33A23 



A10A12 - |17,|2 -A*oA23 - |l^c 



(A.l) 
(A.2) 

(A.3) 



0'33 



A^O (A10A12 - |f}e|2)(A3oA32 - If^cP) 



The trigger susceptibilities are obtained by interchanging the indices {p, 1} o {i,3}. 

For the control field being in the form of a quasi-standing wave fic = fije*'^^^ + ri^e~*'^^^ the medium susceptibilities 
can be then expanded into the Fourier series: 



(1) - ^(1) 



'X.p X' 



p!o + 2 51 Xpln cos{2nkcz), 



and similarly for X^pt2n' ^^'^ explicit form of the Fourier coefficients we rewrite the expressions (jA.ll - 



(1) _ N\dio\^aii AA 



y(3) 
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Xpt 



N\d 



101 



heo 

TjA 



{I + B cos{2kcz)f 



l + Bcos{2kcz) l + Dcos{2kcz) 
if A rj'C* 



l + Bcos(2kcz) 1 + D* cos{2kcz) 
X {S^ + 2ntn- cos(2fccz))] , 

^^10 



TVldioP i^' + 2n+n- cos2kcZ 



heo [1 + B cos{2kcz)Y 



where: A = (^^ - A10A12) \ B = 2n+n-A, S'l = nf 
, C - (52 + A^oA23)"' , D = 2Q+n-C. 



Ai 



' A30A32-A10A 

Let us denote: 



(T33 A12 A23/A13 

" A;oA23+AioAi2 ' 



Ft\,) 



cos(nx)dx 1 

2tT-ttJ 1+^COs(x) yJT~^J? 



1 
1 

1 

1 



dx 



[1 + /iC0s(a:)]2 yr~r^ 

COii{x)dx —fl 



[1 + //cos(a;)]2 



3-4 
^2 



cos(2x)(ia; 2 
[1 + /iC0s(a;)]2 ^2 ' yf— ^ 

cos(3x)(ia; 8 ^ ^ 7; 



3^1 



2t:-^J [1 + //cos(a;)]2 ^3 _ 
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The above integrals have been calculated with the use of the residues method. The Fourier coefficients of the 
susceptibilities read then: 



(1) 

(3) 



N\d,o\^ 
heo 

,Aaf^'\b)-cf^'\d) 



ri'St[AF^'\B)~C*F^''{D 



(1)^ 



(3) 
Xpt,2 



i2ntn- [af[^\b) - c*f[^\d*) 



vf[^\b) 



heo 



- ri[AFl''{B)-CFl'>{D) 



+ iS't [Af[^\b) - C*Fl''{D 
+ ^'ntn-A(F^'\B)+F^'\B))- 

- ri'ntvi-c* (f^^\d*) + f^^\d*) 



(3) 
XptA 



heo 



vFi^\B) 



v[af^'\b)~cf^'\d) 



ri'Si [AF^^\b) ~C*F^'>{D 



Ml 



rfntn:A[F^'\B)+F;,'>{B) 



(3) 



(3) 



N\dio 




heo 




N\dio 




heo 


A^o 


N\dio 




heo 


A^o 



s!F^'\B) + 2ntn~Fi^\B) 



sIf[^\b) + ( f;,^'{b) + f^^'{b) 



SIf!,^\b) + 0+0- ( F[^'{B) + F^iB) 



As shown above, only these coefficients are present in the pulse propagation equations. Note that the corresponding 
formulae given in Ref. |26j included some misprints. 

APPENDIX B 

To evaluate the population redistribution we proceed as follows. We first use equations 2-4 of the stationary version 
of Eqs. (1) to express croo through ajo and ajj: 



O-QO = 5(Op(Tlo) + ^(cril - (T33), 



711 
2 

722 
2 

733 



5(OcCT2o) - ^(cTii + 0-33), 
9(0ta3o) + 2^(733, 



where ^ = fc = 1, 2, 3. Note that here we do not assume any of the populations negligible. Next we express 



11 



(Tjo through the diagonal elements of cr in the lowest-order approximation with respect to the probe and trigger: 

fi*(CTii - CToo) f^*|f^cP(a'22 - O-QO) 



where D = A 
elements: 



where: 



ClO 



|o,|- 



A* 

^10 A 



D 



((722 — CToo) + 



\VLp\^{(Jli - (Too) \%\^\^c?{a22 - CToo) 



A10A12 - |r!,|2 Ai2i?(AioAi2 - If^cP) 



cao 



A30A32 - |r!cP A32i5(A3oA32 - W) 

f^t (0-33 - CToo ) I fic P (0-22 - o-oo ) 



^30 a; 



Combining the two latter steps we obtain a set of linear equations for the diagonal 



CToo = a(CTll — (Too) — ai((T22 — CToo) + C(ctii — (T33), 

CToo = /^(ct33 - CToo) - /?i(ct22 - CToo) + 2^0-33, 

CToo = 7(cti1 - CToo) + <5(ct33 - (Too) + e(CT22 - CToo), 

CToo + CTu + CT22 + CT33 = 1, 



(B.l) 



Pi 

1 

6 




1 



D*iAi,Ai,-\n,\ 
1 



733 

722 

722 



D*{^lo^s2 








i:i(AioAi2 - 




1 




i:i(A3oA32 - 





(B.2) 



722 



1) -\n,\^\np\^^ 



1 



^12 



The solutions of Eqs. (jB.ip are 

1 

M 



CTll 



A32i?2(A3oA32-|r!cP), 

{(1 + a){/3i6 + /3e) + {3 + 2a + I3)e^ + /3i{-a + 6 + 7)^+ 



+ ai[-(5 + 3^ + 2(,5 + 7)^ + ^(l+7 + 0]}, 

CT22 = ^{(-l + «i)/37+[(-l + /3i)'5 + /3(2-ai+e) + 7(-3 + 2ai+/3i)]C+ 

+ 2(3 - ai - A + + <5 + 7 + e)C^ + 

+ « [(-1 + pi)6 + (3 - /3i)e + 2e((5 + e + + /3(1 + e + 2^)]} , 



(B.3) 



CT33 



M 



{-/3i7 + ai(l + /3)(7 - + [(1 + /3)e + /3i(2 + 5 + 7)] ^+ 



+ a[{l + l3)e + pi{l + S + e)]}, 
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while the value of (Too follows from the probability conservation. The denominator M reads: 

M = /3i((5-7) + /3(e-7 + 4ai7) + ai(/3 + 7-(5) + 

+ [2/3 - ,5 + 4e + 4/3e - 87 + 2ai(l - /3 + 87 + (5) + /3i(2 + 47 + 4,5)] C + 

+ 2(3 - 2ai - 2^1 + /3 + 7 + (5 + 2e)e^ + 

+ a [;3 - (5 + e + 4/?e + /3i(l + 4,5 - 20 + 3^ + 2e(/? + 5 + 3e + 0] • 

Using the above equations with constant initial values of probe and trigger fields and a given value of the control field 
we are able to estimate the population distribution. 
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FIG. 2: (color online). The population distribution in the case of trigger present: an - solid red line, (J33 - dashed green line, 
(Too - short-dashed blue line, (J22 - dotted violet line, for a narrow (a) and wide (b) frequency range. The populations croo and 
(722 are negligible and thus not shown in the case (b). 
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FIG. 3: (color online). Probe field's phase shift 7Jp(5i) vs the probe detuning in the presence (solid red) or absence (dashed 
green) of the trigger for balanced population distributions cth = (T33 = 0.5. Phase shift in the presence (dash-dotted blue) or 
absence (dotted orange) of the trigger when instead the population distribution is given by Eqs. (B.3). The frame (a) refers 
to the transmitted probe phase shift {<Pp{5i)) acquired in the running wave configuration, whereas the frame (b) refers to the 
reflected probe phase shift (iPp (<5i)) acquired in the quasi-standing wave configuration. 
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FIG. 4: (color online). Trigger- induced probe phase shifts Atp^ as a function of the probe detuning for transmission (solid red) 
and reflection (dashed green). The population distribution is given by Eqs. (|B.3[1 while all other parameters are as in Fig. |3Jd. 
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FIG. 5: (color online). A comparison of the trigger-induced probe field's phase shifts (a) ipp{6i), (b) (p~{5i), vs the probe 
detuning Si for different control fields: = 6.67 MHz, trigger present - solid red line, trigger absent - dashed green line, 
SlJ — 10 MHz, trigger present - dot-dashed blue line, trigger absent - dotted orange line (cf. Fig. [3}3, lower plot, in which 
n+ = 4 MHz). 
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FIG. 6: (color online). Transmission (a) and reflection (b) coefficients in the presence of a trigger and for different values of the 
control field n+ = 4 MHz (solid red), Q+ = 5.33 MHz (dashed green), fi^ = 6.67 MHz (dotted blue) when the other coupling 
beam is set to fl^ = 2 MHz. 



